Surface-impedance approach solves problems with the thermal Casimir force between 

real metals 
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The surface impedance approach to the description of the thermal Casimir effect in the case of 
real metals is elaborated starting from the free energy of oscillators. The Lifshitz formula expressed 
in terms of the dielectric permittivity depending only on frequency is shown to be inapplicable in the 
frequency region where a real current may arise leading to Joule heating of the metal. The standard 
concept of a fluctuating electromagnetic field on such frequencies meets difficulties when used as a 
model for the zero-point oscillations or thermal photons in the thermal equilibrium inside metals. 
Instead, the surface impedance permits not to consider the electromagnetic oscillations inside the 
metal but taking the realistic material properties into account by means of the effective boundary 
condition. An independent derivation of the Lifshitz-type formulas for the Casimir free energy and 
force between two metal plates is presented within the impedance approach. It is shown that they are 
free of the contradictions with thermodynamics which are specific to the usual Lifshitz formula for 
dielectrics in combination with the Drude model. We demonstrate that in the impedance approach 
the zero-frequency contribution is uniquely fixed by the form of impedance function and does not 
need any of the ad hoc prescriptions intensively discussed in the recent literature. As an example, 
the computations of the Casimir free energy between two gold plates (or the Casimir force acting 
between a plate and a sphere) are performed at different separations and temperatures specific for 
the regions of the anomalous skin effect and infrared optics. The results are in good agreement with 
those obtained by the use of the tabulated optical data for the complex refraction index and plasma 
£^ model. It is argued that the surface impedance approach lays a reliable theoretical framework for 

. the future measurements of the thermal Casimir force. 
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■M ' I. INTRODUCTION 

03 . 

Considerable attention has been focused recently on the Casimir effect [1] which is a rare manifes- 
tation of the zero-point oscillations of the electromagnetic field at macroscopic scales. The Casimir 
force arises as response to the change of the spectrum of zero-point oscillations when material bound- 
^ ■ aries are present. It acts between the boundaries of these bodies and depends on the parameters of 
their materials, their geometry (including surface roughness), and on the temperature (for a detailed 
information see the monographs and reviews [2-6]). Currently, many precision measurements of 
the Casimir force between metal boundaries have been performed [7-14]. Their results were used 
for constraining hypothetical forces predicted by unified gauge theories of fundamental interactions 
[15-17] and in nanotechnological applications [18,19]. 

With respect to the present state of the art, the theoretical description of the Casimir force calls 
for a careful account of all material properties and other relevant factors. Surprisingly, it was found 
that the calculations of the temperature effect on the Casimir force between real metals of finite 
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conductivity run into serious troubles which have been the subject of much controversy [20-33]. 
The key contradiction is on whether the term of the Lifshitz formula [34,35] related to the zero 
Matsubara frequency for the perpendicular polarized modes of an electromagnetic field contributes 
to the physical quantities and, if so, how much would its contribution be. At the moment there are 
five distinct approaches to the resolution of this problem in the recent literature, namely: 

(a) According to the first approach, proposed in Ref. [20] and supported in Refs. [26,27], in the 
case of real metals the term of the Lifshitz formula with zero Matsubara frequency should be 
calculated by using the Drude dielectric function. As a result, the perpendicular polarized 
modes do not contribute to this term. 

(b) In the other approach [28] a special modification of the zero-frequency term of the Lifshitz 
formula, supplemented by the Drude model, was proposed leading to a nonzero contribution of 
the perpendicular polarized modes. This modification was done by analogy with the prescription 
of Ref. [36] for an ideal metal but it does not coincide with it. 

(c) In the framework of the approaches of Refs. [23-25] the modification of the zero-frequency term 
of the Lifshitz formula was made identical to that for an ideal metal [36]. As a consequence, the 
contribution of the perpendicular polarized modes to the zero-frequency term in Refs. [23-25] 
is nonzero and coincides with that for an ideal metal. 

(d) According to Refs. [21,22,29] the contribution of the perpendicular polarized modes with zero 
Matsubara frequency is nonzero and should be calculated by substituting the free electron 
plasma dielectric function into the unmodified Lifshitz formula. 

(e) Finally, according to the approach of Refs. [30,31], the description of the thermal Casimir force 
can be obtained by using the Leontovich surface impedance boundary condition (recently this 
approach was applied also in Ref. [32]). In doing so, the perpendicular polarized modes give a 
nonzero contribution to the zero-frequency term of the Lifshitz formula prescribed by the form 
of the impedance. 

As to the contribution of the modes with a parallel polarization to the zero-frequency term, there is 
consensus between all these approaches that for real metals it is nonzero and the same as for ideal 
metals. Note also that some of the viewpoints varied with the time in the framework of the above 
approaches (a) - (e). For example, in Refs. [23,24] the approach (c) was considered as an universal 
prescription, whereas in Ref. [32] it is restricted by the range of only cryogenic temperatures, and in 
Ref. [25] by the case of sufficiently large separation distances and by the presence of thin covering 
metallic films. The approach (b), proposed in Ref. [28] to resolve the contradictions arising for 
the Drude model in combination with the Lifshitz formula, was considered later in Ref. [31] as 
unnecessary as the Drude model itself turned out to be irrelevant for the description of the thermal 
Casimir force between real metals. It should be particularly emphasized that the approaches (a) and 
(c) were proved to be in contradiction with thermodynamics [31,33] since they violate the Nernst 
heat theorem. A detailed comparison of all the approaches can be found in Refs. [28,31]. 

The present paper aims to work out quite clearly that the surface impedance approach provides 
an answer to all the complicated problems with the retarded Casimir force between real metals 
at both zero and nonzero temperature. We demonstrate that the main reason why the Drude 
model in combination with the Lifshitz theory had failed to describe the thermal Casimir force is 
the inadequacy of the standard concept of a fluctuating electromagnetic field depending only on 
frequency inside a lossy real metal. Rather than to consider fluctuations inside a metal, the surface 
impedance approach suggests that the effective boundary conditions take into consideration in a 
noncontradictory way the involved reflection properties from the surface of a real conductor. In this 



2 



case no additional prescriptions are needed, and the values of the zero-frequency contributions for 
both parallel and perpendicular modes follow immediately from the explicit form of the impedance 
function. 

On this basis, we present for the first time a derivation of the formula for the Casimir free energy 
and force in a configuration of two parallel plates in the surface impedance approach (in Ref. [30] 
it was presented without proof by the use of a prescription changing the integration over continu- 
ous frequencies for the summation over the discrete Matsubara frequencies). The necessity of this 
derivation results from the fundamental role played by the concept of the surface impedance in the 
theory of the Casimir effect between real metals. The relationship between this formula and the 
Lifshitz formula for the free energy is found (the new formula is obtained by exchanging the reflec- 
tion coefficients, which appear in the Lifshitz formula, for the ones derived in the surface impedance 
approach). Our derivation starts from the free energy of the oscillators and suggests also other 
means to derive the usual Lifshitz formula for the thermal Casimir force between dielectrics. The 
obtained formula is applied to compute the Casimir free energy and force at different temperatures 
and separation distances between the test bodies. To do this, one has to use the impedance functions 
describing the regions of infrared optics, anomalous or normal skin effect depending on the value 
of the characteristic frequency giving the main contribution to the Casimir force. It is shown that 
no contradictions with thermodynamics arise, and no artificial prescriptions for the zero-frequency 
term are needed (we also demonstrate that the predictions of large temperature corrections to the 
Casimir force at small separations and cryogenic temperatures, made in Ref. [32], are in error since 
the impedance of the anomalous skin effect was used in [32] outside its range of application). 

The paper is organized as follows. In Sec. II we demonstrate the inadequacy of the concept of a 
fluctuating electromagnetic field inside lossy real metals and remind the basic facts from the theory 
of surface impedance. Sec. Ill is devoted to the derivation of the electromagnetic oscillation spec- 
trum between two parallel plates starting from the impedance boundary condition. In Sec. IV the 
formula for the Casimir free energy is derived in the surface impedance approach. The new simple 
derivation of the usual Lifshitz formula for the thermal Casimir force between dielectrics is also given 
here. Sec. V contains the calculations of the Casimir energy and force at zero temperature using 
the impedance approach. The results are compared with the previously known ones, obtained by 
the use of the usual Lifshitz formula and the tabulated optical data, and found to be in agreement. 
In Sec. VI the computations at nonzero temperature are performed in the impedance approach at 
different separation distances between the test bodies. They demonstrate good agreement in tran- 
sition regions between the different analytical expressions for the impedance. In Sec. VII the reader 
finds conclusions and a discussion of the relationship between the proposed impedance approach to 
the Casimir effect and the theory of the van der Waals forces valid at small separation distances 
between the test bodies. 



II. THE CONCEPT OF A FLUCTUATING ELECTROMAGNETIC FIELD AND THE SURFACE 

IMPEDANCE 

It is well known that the concept of a fluctuating electromagnetic field works well for the description 
of zero-point oscillations within media with a frequency-dependent dielectric permittivity where no 
real electric current does arise. We will now look at a conductor in an external electric field, which 
varies with some frequency u satisfying the conditions 

Z<<J n (w), /<— , (1) 
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where / is the mean free path of a conduction electron, S n (u) = c/y2nauo is the penetration depth 
of the field inside a metal, a is the conductivity, and vp is the Fermi velocity. Eqs. (1) determine 
the domain of the normal skin effect [37]. In this frequency region the external field leads to the 
initiation of a real current of the conduction electrons. 

The normal skin effect is characterized by the volume relaxation described by the temperature- 
dependent relaxation frequency 7(T). As a result, the mean free path of the conduction electrons 
is also temperature-dependent, / = l(T) = vf/^(T), and increases with a decrease of temperature. 
The interaction of the conduction electrons with the elementary excitations of the crystal lattice 
(phonons) leads to the occurrence of electric resistance and heating of the metal. The dielectric 
permittivity of a metal in the domain of the normal skin effect can be modelled by the Drude 
function 



uj [uj + fj{T)\ 

where ui p is the plasma frequency of the free electron plasma model (u p is temperature- independent). 
Remind that the Drude dielectric function (2) was used in the approaches (a) - (c) (see Introduction) 
in combination with the Lifshitz formula to describe the thermal Casimir force between real metals 
for the frequencies both inside and outside of the region (1). This has led to difficulties including 
the violation of Nernst's heat theorem (see Introduction). 

The physical reason for these difficulties becomes quite clear when one observes that the usual 
alternating electric field with frequencies characteristic for the normal skin effect inevitably leads 
to heating of a metal as it penetrates through the skin layer. By contrast, the thermal photons in 
thermal equilibrium with a metal plate or, much less, the virtual photons (giving rise to the Casimir 
effect) can not, under any circumstances, lead to the initiation of a real current and heating of 
the metal (of course, this is strictly prohibited by thermodynamics). Hence the standard concept 
of a fluctuating electromagnetic field penetrating inside a metal described by the Drude dielectric 
function fails to model virtual and thermal photons in the frequency region (1). As a consequence, 
the Lifshitz formula can not be applied in combination with the Drude dielectric function (2) to 
describe the thermal Casimir force even in the domain of the normal skin effect. 

These arguments are supported also by considering the other frequency regions. At higher fre- 
quencies or larger / (lower temperatures) for most of the metals the anomalous skin effect holds, 
which is characterized by the inequalities 

8 a {u) < /, 5 a {u) < — , (3) 

where the skin depth is given by [37] 

and Sf is the total area of the Fermi surface [in fact, within the inequalities (3) "much less" can 
be replaced by "less" thereby preserving Eq. (4) with a good precision]. In the frequency region of 
Eqs. (3) the volume relaxation is not significant but the connection between the electric field and 
current becomes nonlocal. Because of this, a metal can not be described by any dielectric function 
depending only on the frequency. As with the normal skin effect, this leads to the inapplicability 
of the standard concept of a fluctuating field spreading inside a medium described by e(uj) and 
to the impossibility to calculate the Casimir force on this theoretical basis. If one would like to 
preserve the role of a fluctuating field within the domain of the anomalous skin effect, some nonlocal 
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generalization of this concept is required. [Note that for some metals, especially for alloys, instead 
of Eq. (3), the inequalities vp/uj <C / <C 5 hold, specifying the so-called relaxation domain [38,39]; 
here the space dispersion is also essential.] 
On further rise of the frequency, the inequalities hold 

— < S r < /, (5) 
UJ 

where 5 r = c/u p , which determines the domain of infrared optics (note that the condition Tiuj <C ef is 
also supported where Ef = hu p is the Fermi energy). In this domain the volume relaxation does not 
play any role and the space dispersion is absent. Under the conditions (5) metals can be described in 
the framework of the free-electron plasma model with a frequency dependent dielectric permittivity 

e(u) = l~% (6) 

According to the plasma model, the conductivity is pure imaginary and hence there is no any real 
current or heating due to an electric field penetrating into the metal. Because of this, the standard 
concept of a fluctuating electromagnetic field, as a model for the virtual and thermal photons, works 
well in the domain of infrared optics. Remind that the plasma model dielectric permittivity in 
combination with the Lifshitz formula was used to calculate the thermal Casimir force [see approach 
(d) from the Introduction]. This approach did not meet any difficulties or contradictions with 
the basic principles of thermodynamics and has led to physically reasonable results. The domain of 
infrared optics is followed by the domain of ultraviolet frequencies where metals become transparent. 

As is evident from the foregoing, the standard concept of a fluctuating electromagnetic field, 
penetrating into a metal described by the dielectric permittivity depending only on frequency can 
not be used as an adequate model for the virtual and thermal photons of some frequencies. There 
is a frequency region (the domain of the normal skin effect) where this model is in conflict with the 
basic properties of virtual and thermal photons at equilibrium which, among other things, can not 
lead to heating of a metal. In addition, in the region of the anomalous skin effect and relaxation 
domain a metal can not be described by the dielectric permittivity depending only on frequency. 
Because of this, another theoretical basis is preferred to find the thermal Casimir force between real 
metals different from the one used in the case of dielectrics. Here we show that this basis is given 
by the concept of the surface impedance introduced by M. A. Leontovich [35,40]. 

The fundamental difference of the surface impedance approach from the usual approaches is that 
it permits not to consider the electromagnetic fluctuations inside a metal. Instead, the appropriate 
boundary conditions are imposed taking into account the properties of real metal 

E t = Z(u) [B t x n] , (7) 

where Z[uS) is the surface impedance of the conductor, E t and B t are the tangential components 
of electric and magnetic fields, and n is the unit normal vector to the surface (pointing inside the 
metal). The boundary condition (7) can be used to determine the electromagnetic field outside a 
metal. Note, that the impedance Z{uj) and the condition (7) suggest a more universal description 
than the one by means of e. They still hold under the inequalities (3) where a description in terms 
of the dielectric permittivity s(uj) is impossible. For an ideal metal we have Z = and for real 
nonmagnetic metals \Z\ <C 1 holds [40]. 

The calculation of the surface impedance over the whole frequency axis is based on the kinetic 
theory [35] and is rather cumbersome. However, in the domains of the normal and the anomalous 
skin effect, and also for infrared optics, simple asymptotic expressions follow which are of great help 
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to compute the Casimir force between real metals. Thus, in the domain of the normal skin effect, 
given by Eq. (1), the surface impedance is [35] 
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In the domain of the anomalous skin effect, determined by Eq. (3), the impedance depends on the 
shape of the Fermi surface [41]. For a polycrystalline metal, composed of many single crystal grains 
of different orientations, one obtains an approximately spherical Fermi surface and the impedance 
is [41] 

z.m = gil^i^M , ( 9 ) 

where 5 a (u) was defined in Eq. (4). 
In the domain of infrared optics [see Eq. (5)] the impedance is given by [35] 



Now one can impose the impedance boundary condition (7) on the surface of metal plates, find the 
oscillation spectrum in the space between the plates and calculate the Casimir free energy density 
and force without consideration of a fluctuating electromagnetic field inside the metal. This was 
performed at zero temperature in Ref. [42] (see also Ref. [3]). Another approach, being similar in 
spirit, was used at T = in Ref. [43] where the reflection coefficients in the Lifshitz formula were 
expressed in terms of Z(uf). 

Now, the question arises what expression for the impedance (8), (9) or (10) should be used to 
calculate the Casimir effect. To answer this question, it is good to bear in the mind that the main 
contribution to the Casimir free energy and force is given by the frequency region centered around 
the so-called characteristic frequency oo c = c/(2a), where a is the space separation between the two 
bodies, parallel plates for instance. The value of uo c may belong to the frequency region given by 
Eqs. (1), (3) or (5), with the result that the functions (8), (9) or (10), respectively, should be used 
defining the impedance in the domains of the normal or anomalous skin effect and infrared optics. 

By way of example, for most of the metals (Au for instance) at room temperature the application 
region (1) of the normal skin effect with the impedance (8) extends up to the frequencies of order 
10 12 rad/s. The application region (3) of the anomalous skin effect at T = 300 K is very narrow and 
extends up to around (6 — 7) x 10 13 rad/s. It should be stressed, however, that with the decrease of 
temperature the application region of the normal skin effect practically disappears and the anomalous 
skin effect extends to all frequencies being lesser than 10 12 rad/s. The reason is that I increases and 
S n (uj) decreases with the decrease of temperature. As a result, the first inequality in Eq. (1) breaks 
down, whereas the first inequality in Eq. (3) is satisfied also at smaller frequencies. Finally, the 
impedance of the infrared optics (10) is applicable up to the frequencies of order 0.1u p (for Au, for 
instance, u p = 1.37 x 10 16 rad/s). It should be particularly emphasized that the transition frequency 
between the anomalous skin effect and infrared optics does not depend on the temperature because 
all the parameters in the second inequality of Eq. (3) are temperature independent (at temperatures 
much smaller than the Fermi temperature, which is of order 10 5 K). In Sees. V and VI we will discuss 
with more details which impedance function should be used for the calculation of the Casimir force at 
different separation distances between the test bodies (also the transition regions between different 
impedance functions will be considered). 
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III. ELECTROMAGNETIC OSCILLATION SPECTRUM BETWEEN TWO PLATES IN THE SURFACE 

IMPEDANCE APPROACH 



Here the derivation of the photon eigenfrequencies in the framework of the surface impedance 
approach is presented. They are needed to derive the Lifshitz-type formula for the free energy in 
the case of real metals. 

We consider the configuration of two parallel uncharged metal plates, separated by a distance a, 
at a temperature T in thermal equilibrium. Let their nearest boundary planes be described by the 
equations z = ±a/2. We impose the boundary condition (7) on planes z = ±a/2 and determine the 
eigenfrequencies of the electromagnetic field in the free space between the plates. The solutions of 
the Maxwell equations in vacuum can be found in the form 

E a (t, r) = e p (k±, z) exp(ifcj_r_j_ — iut), 

B a (t, r) = g p (k±, z) exp(ik ± r ± - iut). (11) 

Here r— (x,y,z) = (r±,z), a = {p, k±,u}, k± = (ki,k 2 ) is the wave vector in the plane (x, y), u 
is a frequency, and the index p — ||, _L labels the two independent polarization states (|| stands for 
the electric field parallel to the plane formed by k± and z axis, and _L stands for the electric field 
perpendicular to this plane). From the Maxwell equations the oscillatory equations for the functions 
e p , g p follow, 

e p (k±, z) - q 2 (k±, u)e p (k±, z) = 0, 

9 ;(k ± , z) - q 2 (k ± , co)g p (k ± , z) = 0, (12) 

where q 2 (k^_,uj) = q 2 = k 2 — uj 2 /c 2 , the prime denotes the derivative with respect to z, and also the 
first-order equations 

e' Pj3 (k±, z) + ikxtp^k^ z) + ik 2 e Pi2 (k ± , z) = 0, 

g' p , 3 ( k ±i z ) + ihg P ,i( k ±i z ) + ik 2g P ,2(k±, z) = o (13) 

(lower indices 1, 2, 3 after a comma stand for the projections of vectors e p , g p onto the axes x, y, z, 
respectively). 

Substituting Eqs. (11) into the boundary condition (7), using the Maxwell equations to express 
the magnetic field and taking the direction of the normal into account [n=(0, 0, 1) at the plane 
z = a/2, and n=(0, 0, —1) at z — —a/2], we find at the boundaries z = ±a/2, respectively, 



e pA (fc L ,±|)=±^ 
e Pl2 (^>±i)=± ! f 



ih e Pi3 (k ± , ±^ - e' p l (k ± , ±| 
ik 2 e Pi3 (k ± , ±0 - e' p 2 (k ± , ±|) 



(14) 



The same boundary conditions for g p can be obtained also. 

Now, let us consider separately the cases of parallel and perpendicular polarizations beginning 
with the parallel one. Without loss of generality, we temporarily assume that k 2 = 0. In this case 
e \\,2(k±, z) = and the solution of Eq. (12) has the form 

e \\,i(k±, z) = B sinh qz, e\\,3(k±, z) — B cosh qz, (15) 

where A and B do not depend on z. From Eq. (13) it follows Aq + ik x B = 0. For the sake of 
convenience, we choose 
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(16) 



Substituting Eqs. (15) and (16) into the impedance boundary condition (14), one obtains the dis- 
persion equation for the spectrum of the electromagnetic oscillations between plates: 

A^V, k ± ) ee e--/ 2 (sinh ^ - ^ cosh = 0. (17) 

Eqs. (12), (13) also have the solutions 

e\\ i i(k±,z) = e~ aq l 2 coshgz, e\\,2(k±, z) = 0, 

e p (^^) = -ye- a< '/ 2 sinh^. (18) 

After substitution of Eq. (18) into Eq. (14) a further dispersion equation for the modes with parallel 
polarization is obtained: 

Af(u, k ± ) = e--/ 2 (cosh f-^ sinh |) = 0. (19) 

It is obvious that in Eqs. (17) and (19) = (fci, A^) can be now considered as arbitrary. 

Exactly the same procedure is applicable to the case of the perpendicular polarization. Once again, 
assuming temporarily k 2 = 0, we obtain the solutions of Eqs. (12) and (13) in the form 

e±,i(k ± ,z) = e ±t3 (k ± ,z) = (20) 

and 

e±, 2 (k±,z) = e _a<?/2 sinh or e ±j2 (k±, z) = e~ aq/2 coshqz. (21) 

Substituting these solutions into Eq. (14), we arrive at two dispersion equations for the determination 
of the electromagnetic eigenfrequencies with perpendicular polarization 

A«( W , k ± ) == e--/ 2 fsinh a A + ^ cosh a A) = 0, 
V 2 to 2 J 

A? ( w , k ± ) ee e~-/ 2 (cosh I + ^ sinh f ) = 0. (22) 

Let us denote the solutions of the transcendental equations (17), (19) by uj\ n , and the solutions 

of the transcendental equations (22) by oj^ ±n . Multiplying Eqs. (17), (19), we can finally find oj\ ±n 
from the equation 

A||(w, k ± ) = A<[V k ± ) Af(u, k ± ) (23) 
= ^ e ^ aq (l - V 2 ) (sinhag - 2 coshagj = 0, 

where 77 = r)(u) = Zuj/{cq). 

In perfect analogy to this, by multiplication of Eqs. (22), one can find the equation for the deter- 
mination of OJ^ \ 
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A ± (cu, k ± ) = A? (a;, k ± ) A ( j°(w, k ± ) 

= - e" a<? (l - k 2 J f sinhag + coshagj = 0, 



(24) 



where k = k(u) = Zcq/uo. 

Note that we have obtained the conditions for the determination of the electromagnetic oscillation 
spectrum by the use of equations for e p . Exactly the same spectrum is obtained if the equations for 
g p are used. 

IV. CASIMIR FREE ENERGY IN THE SURFACE IMPEDANCE APPROACH 



Now we are in a position to present a rigorous derivation of Lifshitz-type formulas for the Casimir 
free energy and force for the configuration of two plates at a temperature T in thermal equilibrium in 
the surface impedance approach. As shown below, these formulas are well adapted for the calculation 
of the Casimir effect between real metals and are not subject to the disadvantages of the approaches 
(a) - (d) discussed in the Introduction. 

First we consider the case of real eigenfrequencies u\ n , 0J^ ±n (this is fulfilled for the pure imaginary 
impedance). The total free energy of the electromagnetic oscillations is given by the sum of the free 
energies of separate oscillators over all possible values of their quantum numbers, 



— - + k B T\n M - e ^ T 



where k B is the Boltzmann constant. Identically, Eq. (25) can be rewritten as 



T = k B T s £\n (2 sinh 

a \ 



2k B T 



(25) 



(26) 



It is clear that at T — > 0, the value of T from Eqs. (25), (26) coincides with the sum of the zero-point 
energies which is the traditional starting point in theoretical investigations of the Casimir effect at 
zero temperature. 

Applying this to the electromagnetic oscillations between metal plates, where a = {p,k±,n}, we 
obtain 



Jo 2n „ 



ln(2sinh5^] + In ( 2 sinh % ^ 
2k B T { 2k B T , 



(27) 



According to the calculations of Sec. Ill, the eigenfrequencies of the electromagnetic field between 
plates with parallel and perpendicular polarizations are determined by Eqs. (23) and (24), respec- 
tively. 

The expression in the right-hand side of Eq. (27) is evidently infinite. Before performing a renor- 
malization, let us equivalently represent the sum over the eigenfrequencies uj\^ n by the use of the 
argument theorem like it is usually done in the derivation of the Lifshitz formula at zero temperature 
by the method of surface modes [6,44,45]. Then Eq. (27) can be rewritten as 



J==knT 



I 

Jo 



00 k\ dk\ 1 



2ti 2m JCx 



In I 2 sinh —7—7= I d 
V 2k B T J 



lnA||(w, k±) +lnA ± (u,k 1 _) 



(28) 



Here, the closed contour C 1 is bypassed counterclockwise. It consists of two arcs, one having an 
infinitely small radius e and the other one an infinitely large radius R, and two straight lines L 1 , L 2 
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inclined at the angles ±45 degrees to the real axis (see Fig. la). The quantities A\\ t ±(u, k±), having 
their roots at the photon eigenfrequencies, are defined in Eqs. (23), (24). Note that, unlike the usual 
derivation of the Lifshitz formula at nonzero temperature [2] , the function under the integral in (28) 
has branch points rather than poles at the imaginary frequencies u\ = where 



2rck B Tl 



I = 0, ±1, ±2, 



(29) 



are the Matsubara frequencies. The contour C\ in Fig. la is chosen so as to avoid all these branch 
points and to enclose all the photon eigenfrequencies. 
The integral in Eq. (28) can be calculated as follows: 



1,-L 



= <i> In 2 sinh 



2m Jci 
1 



hu 

2~k B T / 



din A|| j j_(a;, k±) 



(30) 



2m 



/ In (2 sinh ) din Aim (u, k ± ) + [ In ( 2sinh ) din Aim (u, k± 

Jl 2 \ 2k B T) ll,±V ' ±J Jc R V 2k BTJ 



/ 



+ / In 2 sinh 



hu 

2k B T / 



dhiA\\ t ±(u, k± 



In 2 sinh 



2kBT, 



din Aim (a;, k±) 



The integral along the arc of infinitely large radius C R vanishes which follows from Eqs. (23), (24) 
under the natural conditions 



lim Z{uj) = const, 



, dZ(uj) 
lim — = 0. 



du 

Integrating by parts in the right-hand side of Eq. (30), one obtains 



(31) 



2iri 



In \ '1 si'ih 1 In A|| ; j_ (a,?, A;_j_) 



+ In ^2 sinh 
+ In 2 sinh 



huj 
2k B T / 

hu 

2k B T / 



2k B T / 
lnA l± (u,k ± ) 

In A|| ! _ L (o;, k±) 



h 



h 



B k B T J Li 

- i£ h 



k B T Jl 2 

hu 

coth 



k B T Jc, 



la 



coth 



2k B T 

hu 



2k B T 



2k B T 
In A|| )J _(a;, k±)du 

In Aim (a;, k±)du 



(32) 



where the contours L 12 and the points A, B are shown in Fig. la. It is evident that all terms, 
besides the integrals, cancel each other or are equal to zero (at the points A, B). The integral along 
the Li can be calculated by the application of the Cauchy theorem to the closed contour C 2 (see 
Fig. lb) inside of which the function under consideration is analytic 



/ coth— — — In Aii ±(u, k± )du = — [ coth— — — In Am ±(u, k±)du. 
J Ll 2k B T v ; Jioo 2k B T v ; 



(33) 



Here it is taken into account that the integral along Cr vanishes. The path (ioo; is) contains 
semicircles of radius e about the singular points (poles) of the function coth(hu /2k B T) . The 
analogous formula for the integral along the line L 2 is 



L 



hu 

coth — — — 
l 2 2k B T 



In 



hu 

coth — — — In A|| i± (a;, k ± )du. 

-ie ZK B 1 



(34) 



Substituting Eqs. (33), (34) into Eq. (32), one arrives at 
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h± = 



coth — — — In A|| 5 j_(cj, k±)dou. 



2irikBT Jioo 2k B T 
The integration in Eq. (35) involving poles at the points i& leads to 



in 



2nk B T 

00 

n ^2 res 

i=— 00 



2k bT 



inA|, ± (e,fcL)de 



(35) 



(36) 



coth In A|| i± (w, k±);i£i 



where the functions A\\ ; ±({;,k±) are obtained from A\\ i ±(u,k±) by the substitution u = i£. In the 
case of real eigenfrequencies, which is under consideration now, Aim are even functions of u (and £). 
As a consequence, the seemingly pure imaginary integral in the right-hand side of Eq. (36) vanishes. 
After the calculation of the residues, and using the evenness of the functions A^±(u, fcjj, the result 
is 



V = E' lnA iuteA), 



(37) 



1=0 



where the prime on the summation sign means that the term for / = has to be multiplied by 1/2. 

Substituting the values (37) of the integrals (30) into Eq. (28), we find the equivalent but more 
simple expression for the Casimir free energy 



T=^-\ k ± dk ± Y, hiA,|te,A; ± )+lnA ± te,A; ± ) 



(38) 



Expression (38) is still infinite. To remove the divergences, we subtract from the right-hand side 
of Eq. (38) the free energy in the case of infinitely separated interacting bodies (a — > 00). Then the 
physical, renormalized, free energy vanishes for infinitely remote plates. From Eqs. (23), (24) after 
the substitution u — > i£i in the limit a — > 00 it follows 



a?(&*0 = j(i + «?) (i + T^)- 

The renormalization prescription is equivalent to the change of Ay ,j_(£z, k±) in Eq. (38) for 



(39) 



A?,(6A) 



(40) 



where the quantities r\\ t ±(£t, k±) have the meaning of the reflection coefficients and are given by 



i(&A) = 



I -Vi 

l + rn. 

1 - Kl 
l + Ki 



' cqi - Zijx 
cqi + Ziii 



1CQ1X 



6 + Z t cqi I 



(41) 



Here Z\ = Z(i£i) and qf = k]_ + The reflection coefficients (41) are in accordance with Ref. [40] 

where the reflection of a plane electromagnetic wave incident from vacuum onto the plane surface of 
the metal was described in terms of the surface impedance. 
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In such a manner the final renormalized expression for the Casimir free energy in the surface 
impedance approach is given by 

Fn = ^~ k ± dk ± £ 1 - rjj(& ^)e" 2 ^ + In 1 - rife, k^e' 2 ^] } , (42) 

where the reflection coefficients are given by Eq. (41). 
The Casimir force, acting between plates, is obtained from Eq. (42) 



dT R k R T 



da 

x 



B± f°° I 

— / k ± dk ± y] q t 

7T Jo l=Q 



(43) 



rfteA)^ - l]" 1 + [r: 2 (6A)e 2a * - l] -1 } 



The above derivation was performed under the assumption that the photon eigenfrequencies are 
real. This is, however, not the case for arbitrary complex impedance. If the photon eigenfrequencies 
are complex, the free energy is not given by Eq. (26) (which is already clear from the complexity 
of the right-hand side of this equation). For arbitrary complex impedance the correct expression 
for the free energy should be determined from the solution of an auxiliary electrodynamic problem 
[46]. It turns out that the Casimir free energy and force are the functionals of the impedance even 
when the impedance has a nonzero real part taking absorption into account. The solution of the 
auxiliary electrodynamic problem leads to conclusion [46] that the correct free energy is obtained 
from Eqs. (38)-(42) by analytic continuation to arbitrary complex impedance, i.e., to arbitrary 
oscillation spectra. The qualitative reason for the validity of this statement is that the free energy 
depends only on the behavior of Z{uj) at the imaginary frequency axis where Z{u) is always real [see, 
e.g., Eqs. (8) - (10)]. Note that in the case of complex eigenfrequencies exactly Eqs. (42), (43) should 
be used written in terms of summations from zero to infinity. Although for real eigenfrequencies 
the summations over / from — oo to oo can be equivalently used it is not so for complex u a as the 
dispersion functions A\\,± cease to be even any more [46]. 

It is necessary to stress that the above derivation of the free energy in the impedance approach can 
be simply modified in order to present the new derivation of the usual Lifshitz formula describing the 
thermal Casimir force between dielectrics. In fact, nothing should be changed in the presentation of 
this section except for the explicit expressions of the dispersion functions A\\,± in Eq. (38) and thus of 
the reflection coefficients r\\ t ± in Eqs. (42) and (43). The dispersion functions should be determined 
not according to Sec. Ill but from the consideration of a fluctuating electromagnetic field both inside 
and outside of the dielectric plates with the usual boundary conditions at the interfaces [2,6,44,45]. 
As a result, the Lifshitz reflection coefficients take the form 

^.w-feav ^-(f^)". <44> 

where S\ = e(u) is the dielectric permittivity of the plate material, and kf = k]_ +e^f/c 2 . 

Then the Lifshitz expressions for the Casimir free energy and force between dielectrics are given 
by Eqs. (42), (43) where the substitution ry, r± — > rii £,, r±^ is made. In such a manner, we have 
performed also a new derivation of the usual Lifshitz formula between dielectric plates starting from 
the free energy of an oscillator. Conversely, the free energy and force in the framework of the 
impedance approach are obtained from the Lifshitz formula if the Fresnel-type reflection coefficients 
r \\,L, r ±,L are changed for those obtained by the use of the impedance boundary condition. 

It should be stressed, however, that the reflection coefficients (44) differ essentially from the 
impedance coefficients (41). To take an example, it is not possible to obtain the coefficients (41) 
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from (44) even if both descriptions in terms of e(uj) and Z(u) are applicable and the impedance is 
expressed in terms of the dielectric permittivity by means of the relation Z(u) = l/yje(u) (which 
holds, e.g., in the region of infrared optics). This underlines the fundamental role of the impedance 
boundary condition as an alternative to the consideration of a fluctuating field inside a medium 
described by e(oo) in the case of real metals. 

We conclude this section by remarking that the obtained expression (42) for the free energy gives 
the possibility also to find the thermal Casimir force in configuration of a sphere (spherical lens) 
above a plate made of real metals in the surface impedance approach 

F(a) = 2irRF R (a), (45) 

where R is the sphere radius. The approximate expression (45) is obtained by the application of the 
proximity force theorem [6] and has an accuracy around a fraction of a percent for configurations 
used in precision experiments on the measurement of the Casimir force [7-12,14,18]. Thus, the 
impedance approach provides the theoretical basis for the measurements of the thermal Casimir 
force between real metals to be performed in near future. 



V. CALCULATION OF THE CASIMIR ENERGY IN THE SURFACE IMPEDANCE APPROACH 



First, we apply the obtained general formulas at zero temperature. In this case Eq. (42) for the 
free energy transforms to the double integral representing the Casimir energy between plates [or, 
according to Eq. (45), the Casimir force acting between a sphere and a plate] 



ii P OO P OO 

E(a) = ^J Q k ± dk L J q di {in [l - rj(£, k^e^] + In [l - ri(£, k ± )e~ 2a «} } , 



(46) 



where the reflection coefficients in terms of the impedance are given by Eq. (41) with the substitution 

(47) 



Let us calculate the quantity (46) obtained in the impedance approach and compare the results 
with the available data found by the traditional computations using the Lifshitz formula. For the 
purpose of numerical computations, it is convenient to rearrange Eq. (46) to the form [30] 



E(a) 



he 



P OO P OO \ 

J d(J ydy 21n(l-e^)+ln 



1 + 



ey 



+ In 



1 + 



ey 



32ir 2 a 3 Jo J( 
where the dimensionless variables (, y are defined as 

u c c 

and the quantities XH'- L (^,|/) are given by 

II(( ' !)= (FW' xHCv) 

The first contribution in the right-hand side of Eq. (48) describes the case of an ideal metal 



V = 2qa, 



Kyz 



720a 3 ' 



(48) 



(49) 



(50) 



(51) 
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the others are the corrections due to the finite conductivity. 

As was stressed in Sec. II, with the decrease of temperature the range of application of the normal 
skin effect (1) reduces to zero, and at T = only the anomalous skin effect and infrared optics occur 
with the frequency regions given by Eqs. (3), (5), respectively. The transition frequency Q between 
the two effects can be obtained from the equations 

= 77 = 8r = —, (52) 

where, according to Eq. (4), S a (Q) = Ca/fi 1 / 3 . All computations given below are performed for Au 
withcjp = 1.37 x 10 16 rad/s [47] and-u^ = 1.4xl0 6 m/s (see, e.g., [48]). Then from Eq. (52) we obtain 
the values of both C a = 8.8 x mrad 1//3 /s 1 / 3 and Q = 6.36 x 10 13 rad/s. If to consider Q as the 
characteristic frequency giving the main contribution to the Casimir effect (Q = u c = c/2a tr ), the 
transition separation distance between the two effects turns out to be equal to a tr = 2.36 /an. Then 
it follows that at distances A p < a < a (r = 2.36 /im the impedance of the infrared optics determines 
the value of the Casimir energy and force, whereas at a > a tr = 2.36 /zm the impedance of the 
anomalous skin effect is applicable (A p =137nm is the plasma wavelength for Au). Direct calculations 
by Eqs. (48), (50) show that the main contribution to the Casimir energy is given by the narrow 
frequency interval around the characteristic frequency u c . Thus, the interval (0.1cj c , 10o> c ) contributes 
94% of the total energy in the wide separation region. What is even more important, the remainder 
does not depend on the form of the impedance function outside of the interval (0.1u; c , 10u c ), to within 
the error of about 0.5%. From this it follows that at each separation distance between the plates 
one should, first, determine the characteristic frequency uj c and, second, fix the proper impedance 
function. Thereafter the chosen impedance function can be used at all frequencies when performing 
the integration in Eq. (48). At zero temperature this prescription is optional. At T ^ 0, however, it 
takes on great significance (see Sec. VI). 

In Fig. 2 the correction factor to the Casimir energy E{a)/E^\a) is plotted which is computed 
by Eqs. (48), (50) and (51) as a function of the separation distance. The solid line is obtained 
with the impedance of the infrared optics (10), and the dotted line with the impedance of the 
anomalous skin effect (9). Both lines are plotted at all separations a > X p to make sure that each 
impedance function is applicable within its own frequency region and to follow their applicability 
at the transition separations around a tr . It must be emphasized that the solid line coincides with 
the correction factor to the Casimir energy computed on the basis of the usual Lifshitz formula in 
combination with the dielectric function of the plasma model (this was demonstrated in detail in 
Ref. [30]). Thus, both the impedance approach and the Lifshitz formula combined with the plasma 
model, lead to one and the same result for the Casimir energy at separations a > \ p . 

As is seen from Fig. 2, at A p < a <C 2.36 jum, the pointed line computed with the impedance of the 
anomalous skin effect (which is inapplicable in this region) significantly underestimates the correction 
factor due to the finite conductivity. For example, at a = 0.15 /xm the values of the correction factors, 
given by the solid and dotted lines, are 0.623 and 0.851, respectively, i.e., the error introduced by the 
use of the impedance of the anomalous skin effect is almost 37%. At a separation a = 0.5 /im this 
error is more than 9%, and decreases with increasing separation. Notice that the computations on 
the basis of the usual Lifshitz formula and optical tabulated data for the complex refraction index 
[which are used to obtain e(i£) through the dispersion relation] also practically coincide with those 
given by the surface impedance in the region of the infrared optics (solid line in Fig. 2). Thus, at 
the separations of 0.2 /im, 0.5 /im and 3 /im the correction factor obtained by the tabulated data 
and Lifshitz formula is equal to 0.69, 0.85, and 0.97, respectively [45,47], whereas in the impedance 
approach it takes the values 0.689, 0.849, and 0.972. 

At larger separations (a ^> a tr = 2.36 /im) the impedance function of the anomalous skin effect 
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should be used to compute the Casimir energy (dotted line in Fig. 2). As is seen from that figure, at 
these separations the impedance of the infrared optics overestimates the role of the finite conductivity 
corrections to the Casimir energy. This overestimation is, however, to within a fraction of a percent. 
In the transition region a = 2—2.5 /xm the results given by both impedance functions are in agreement 
bringing the discrepancies of about 1% only. This leads us to the conclusion that at zero temperature 
both impedance functions work well in their respective application regions. In the transition region 
each of them can be applied and the results are in agreement within an error of 1%. It is seen also 
that "much less" or "much larger" in the above inequalities in fact mean two or three times less 
(larger). 

If to speak about the region of infrared optics, the Lifshitz formula in combination with the plasma 
model or optical tabulated data for the complex refraction index leads to the same results as the 
impedance approach. It gives rather good results even in the region of the anomalous skin effect, 
where, strictly speaking, the description in terms of e is not applicable (see Sec. II). The feasibility 
of the Lifshitz formula is explained by the fact that at zero temperature the normal skin effect is 
practically absent and the problems connected with the heating of a metal due to the real electric 
current are not relevant. As a result, both the impedance approach and the usual Lifshitz formula 
are applicable. At nonzero temperature, however, the surface impedance approach acquires a new 
meaning and solves the problems formulated in Introduction (see the next section). 



VI. CALCULATION OF THE CASIMIR FREE ENERGY IN THE SURFACE IMPEDANCE APPROACH 



Here we calculate the Casimir free energy for the configuration of two parallel plates made of Au at a 
temperature T at thermal equilibrium. The starting point is Eq. (42) where the reflection coefficients 
are expressed in terms of the surface impedance by Eq. (41). Introducing the dimensionless variables 
by analogy with Eq. (49), we transform Eq. (42) to a form convenient for numerical computations: 



+ ln 



1 + 



+ In 



1 + 



(53) 



where X", X 1 - are given by Eq. (50) with the change ( — > Q, Z — > Z t — Z(ic(i/2a). Notice that the 
first contribution in the right-hand side of Eq. (53) presents the Casimir free energy for the ideal 
metal [6] 



(54) 



i=i 



+ 7T 



T 



T, 



eff, 



^sinh 2 [nl T 




Icoth^ 



T, 



eff, 



where E^(a) is defined in Eq. (51) and the effective temperature ksTeff — hw c — he/ (2a). 

First of all, let us demonstrate that in the impedance approach there is no problem with the 
contribution of the zero Matsubara frequency which was the subject of much recent controversy 
(see Introduction). We start from the lowest characteristic frequencies where the impedance of the 
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normal skin effect, given by Eq. (8), is applicable. Substituting it into Eq. (41) and putting £j = 0, 
one obtains 



rj(0A)=ri(0A) = l, (55) 

i.e., the same result as for an ideal metal. Namely, this becomes clear since the quantities X^ 1 , 
defined in Eq. (50) and given for the impedance of the normal skin effect as 

X"(0 j| /)=X- L (0,y)=0, (56) 

when inserted into Eq. (53) obviously lead to the same zero-frequency contribution as it holds for 
an ideal metal. It should be stressed that all functions r?\_ L (£,k±) and XH'- L (^, y) are continuous 
functions of two variables including the point (0,0). Thus, the case of an ideal metal is achieved as 
a limiting case of a real metal with increase of the conductivity when the real metal is described in 
the framework of the impedance approach. Remind that this is not the case when the real metal 
is described by the Drude dielectric function (2) and the Lifshitz formula for dielectrics is used to 
calculate the Casimir free energy [approach (a) from Introduction]. In fact, if doing so it follows 
from Eq. (44) 

rf L (0A) = l, ri >L (0A) = 0, (57) 

and there is a break of continuity between the properties of real metals and of ideal metal [28]. 

At higher characteristic frequencies the anomalous skin effect holds with an impedance function 
as from Eqs. (4), (9). If to extend this function to all frequencies (to zero Matsubara frequency in 
that case) we ensure that both the Eqs. (55) and (56) are valid once again. As a result, in both 
regions of the normal and the anomalous skin effect the thermal corrections to the Casimir free 
energy and force for real metals are very close to those for an ideal metal. As one would expect, 
at large separations (characteristic for the anomalous and especially for the normal skin effect) all 
metals behave like an ideal one [this is, however, not the case in the framework of the approach (a)]. 

If the characteristic frequencies increase further, the infrared optics with an impedance function 
of Eq. (10) takes place. The extension of it to zero Matsubara frequency leads to 

i«UJ = l, ri (0,y=(^|) 2 . (58) 

In this case there occurs a dependence of the perpendicular reflection coefficient at zero frequency 
on the properties of the real metal through the value of the plasma frequency. This is reasonable, 
because the real properties of a metal are most pronounced at small separations characteristic for 
the infrared optics. In the limit uo p — > oo the result for an ideal metal is reproduced from Eq. (58). 

Before performing the computations, it must be emphasized that the surface impedance approach 
is in perfect agreement with thermodynamics. In the impedance approach the entropy, defined as 

S (a ,T) = -^l>, (59) 

is positive and equal to zero at zero temperature in accordance with the Nernst heat theorem [remind 
that this is not the case in the approaches (a) and (c)]. The validity of the Nernst heat theorem 
in the impedance approach can be demonstrated in the regions of both the infrared optics and the 
anomalous skin effect (as noted above, the region of the normal skin effect dies out with decreasing 
temperature). According to the results of Ref. [30], in the region of infrared optics the Lifshitz 
formula combined with the plasma model leads to exactly the same perturbation results for the 
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Casimir free energy and force as the impedance approach. At T T e ff the free energy is given by 
[31] 



T R {a,T) = E{a) - 



ft<(3) 
167ra 3 




(60) 



where E(a) is the Casimir energy at T = defined in Eq. (46). After the substitution in Eq. (59), 
this leads to the simple expression for the Casimir entropy 



S(a,T) = 



3k B ((3) ( T 
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eff, 
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135C(3)T e// 



(61) 



which is positive and equal to zero at zero temperature. 

The impedance approach in the region of the anomalous skin effect was used in the recent Ref. [32]. 
The asymptotic of the entropy at very low temperatures, obtained in Ref. [32] , demonstrates that it 
is positive and has zero value at zero temperature in accordance with the requirements of thermo- 
dynamics. 

By the way of example, here we perform the numerical computations of the relative thermal 
correction to the Casimir free energy defined as [J r R(a, T) — E(a)]/E(a). This quantity has also 
the meaning of the relative thermal correction to the Casimir force in the configuration of a sphere 
(spherical lens) above a plate used in precision experiments on the Casimir effect. If the characteristic 
frequency u c belongs to the region of the normal skin effect, the results practically coincide with 
those obtained for an ideal metal [30], and the free energy is given by Eq. (54). If the characteristic 
frequency belongs to the regions of the anomalous skin effect or infrared optics, the computational 
results for the relative thermal correction are obtained by Eqs. (48), (53) and presented in Fig. 3. 
The solid lines are computed with the impedance of the infrared optics (10), and the dotted lines 
with the impedance of the anomalous skin effect (9). All computations are performed for Au at two 
temperatures T = 300 K and T = 70 K with numerical parameters as listed in Sec. V. Both pairs of 
lines are plotted at all separations a > \ p for a better visualization of the application range of each 
impedance function. 

Remind that at separations between the plates X p < a < a tr , where a tr = 2.36 /im does not 
depend on the temperature, the impedance of the infrared optics is applicable, and at separations 
a > a tr the impedance of the anomalous skin effect should be used. It is seen from Fig. 3 that 
at small separations the use of the impedance function of the anomalous skin effect significantly 
overestimates the value of the thermal correction. Thus, at a = 0.15 /xm the values of the relative 
thermal corrections given by the dotted and solid lines are 1.55 x 10~ 2 and 1.82 x 10~ 4 , respectively, 
at T = 300 K, and 4.85 x 10" 3 and 2.76 x 10~ 6 , respectively, at T = 70 K. What this means is the 
thermal correction, predicted by the impedance of the anomalous skin effect in the region of the 
infrared optics, where this impedance is not applicable, is in 85 times greater at T = 300 K and in 
1757 times greater at T = 70 K than the correct values. 

At separations a > a tr = 2.36 /zm the dotted lines present the correct dependence of the thermal 
correction on the separation distance. The difference between the free energies computed by the use 
of two impedance functions is, however, to within a fraction of a percent. In the transition region 
the results, given by the impedance function of the infrared optics and anomalous skin effect, are in 
agreement with a sufficient accuracy. For example, at a = 2.5 /xm the ratio of the relative thermal 
corrections obtained by the use of different impedance functions is 1.05 at T = 300 K and 2.19 at 
T = 70 K. As a result, the discrepancies in the values of the free energy are about 1.2% (T = 300 K) 
and 0.7% (T = 70 K). 
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Our results for the thermal correction to the Casimir free energy are in disagreement with the 
conclusion of Ref. [32] about the existence of large thermal corrections at low temperature made 
in the framework of the impedance approach. As correctly argued in Ref. [32], the description of 
metals with the impedance in the region of the anomalous skin effect is more appropriate than with 
the dielectric permittivity. However, the conclusion about the existence of large thermal corrections 
at separations 100 nm-500 nm at T < 70 K made in that paper is in error. To obtain this conclusion, 
the impedance function of the anomalous skin effect was applied in Ref. [32] at separations much 
less than a tr = 2.36 /im, i.e., in the separation range of the infrared optics. This was explained by 
the fact that at temperatures T < 70 K the inequality I ^> 5 r holds which, from the standpoint of 
Ref. [32], guarantees the applicability of the impedance of the anomalous skin effect. In actual truth, 
this inequality is not sufficient. In fact, one additional inequality, S a (u) <C vp/cu, must be fulfilled in 
order that the anomalous skin effect holds [see Eq. (3) and Ref. [37]]. Because of this, the frequency 

[see definition in Eq. (52)], considered in Ref. [32] as the characteristic frequency of the anomalous 
skin effect, is actually the transition frequency to the region of infrared optics. As a result, all 
computations performed in Ref. [32] correspond to the dotted line at T = 70 K of our Fig. 3 at 
separations a < a tr = 2.36 /im. According to our computations, in this separation range the dotted 
line at T = 70 K overestimates the value of the thermal correction by a factor of 2000, whereas the 
correct results are given by the solid lines obtained by the use of the impedance function of infrared 
optics. Remind that the characteristic frequencies, corresponding to the separations 100 nm-500 nm, 
fall into the interval u c = (0.3 — 1.5) x 10 15 rad/sS> fl, i.e., belong to the region of the infrared optics 
(see Sec. II). 

At the end of this section, we would like to stress that in the sums, like in Eqs. (42), (43) and (53), 
the form of the impedance at the characteristic frequencies must be substituted and extended to all 
other frequencies. At zero temperature, as was shown in Sec. V, the frequency region [0, lu c , 10u; c ], 
where the characteristic frequency is oj c = c/(2a), gives most of the contributions to the result. 
Calculations show that at nonzero temperature the Matsubara frequencies from £ to £at ~ 10co» c 
give the dominant contribution. For example, at a = 0.15 /xm (uj c = 10 15 rad/s), T = 300 K the first 
41 Matsubara frequencies determine the total result. Here £i = 2.5 x 10 14 rad/s and £40 = 10 16 rad/s. 
All nonzero Matsubara frequencies belong to the region of infrared optics. With a decrease of the 
temperature some of the Matsubara frequencies may fall within the frequency region of the anomalous 
skin effect (at T = 70 K, for instance, the first Matsubara frequency £1 = 5.75 x 10 13 rad/s< f2 = 
6.36 x 10 13 rad/s for Au) . At small separations, however, the differences in the contributions of several 
first Matsubara frequencies computed by the use of different impedance functions are negligible. At 
T = any extension of the impedance function outside of the above interval leads to approximately 
one and the same value of the Casimir energy (in the integral one point £ = is of no significance). 
At T ^ 0, however, the contribution of the zero Matsubara frequency £0 — becomes dominant 
at large separations (high temperatures), and at room temperature, for instance, it determines the 
total value of the free energy at a > 5 /im. 

The basic challenge is whether the actual reflection properties of plate materials at very low, 
quasi-static, frequencies are responsible for the Casimir force in the high-temperature limit. The 
point is that materials are at hand (e.g., indium tin oxide) which are good conductors at quasi-static 
frequencies but transparent to visible and near infrared light. If to consider a pair of plates made of 
indium tin oxide (ITO) at a separation a = 5 /im, and the other pair of plates at the same separation 
made of Au, one runs into difficulties. If the actual low-frequency reflection properties should be 
substituted into the zero-frequency term, the impedance of the normal skin effect from Eq. (8) must 
be used. As a result, the thermal Casimir force at a = 5 /im will be equal for both pairs of plates 
(and practically the same as for an ideal metal). This is in contradiction with the physical intuition 
as around the characteristic frequency u c = 3 x 10 13 rad/s (computed at the separation 5 /im) ITO 
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is a poor reflector. A better physical result would be obtained if one extends the characteristic 
impedance at u c (of the anomalous skin effect for Au and of the infrared optics for ITO) to zero 
Matsubara frequency. If this is done, the zero-frequency term for Au plates will be the same as for 
an ideal metal in accordance with Eq. (55). For ITO plates the zero-frequency term will contain the 
value of u>p TO according to Eq. (58). Taking into account the large value of the penetration depth for 
ITO, the magnitude of the Casimir force between the ITO plates will be less than between the plates 
made of Au, as the intuition suggests. In fact, the question on whether the values of the Casimir 
force for the two above pairs of plates at separation 5/im are equal or different, can be answered 
experimentally using the measurement scheme suggested recently in Ref. [49]. We expect that the 
experimental result will be in accordance with the suggestion of the physical intuition (note that 
this example with two pairs of plates was used with another aim in Ref. [31]). 

VII. CONCLUSIONS AND DISCUSSION 

In the foregoing we have presented the surface impedance approach to the theory of the Casimir 
effect at both zero and nonzero temperature. In the impedance approach the effective boundary 
condition is imposed taking the real properties of the metal into account. Previously this approach 
was considered as nothing more than a useful approximation to the more complete Lifshitz theory 
using the concept of a fluctuating electromagnetic field both outside and inside the boundary of the 
bodies. Our conclusion is that the standard concept of a fluctuating field inside a metal, described 
by the dielectric permittivity depending only on frequency, in the region where a real current may 
arise, can not serve as an adequate model for the zero-point oscillations and thermal photons. It 
follows from the fact that the vacuum oscillations and the thermal photons in equilibrium under 
no circumstances can lead to a heating of the metal. If this fact is overlooked, contradictions with 
the thermodynamics arise when one substitutes into the Lifshitz formula for the Casimir free energy 
and force the Drude dielectric function which takes into account the volume relaxation and, conse- 
quently, the Joule heating. This situation reflects the nontrivial character of quantum fluctuations 
in nonhomogeneous case involving both vacuum and real metals, containing conduction electrons, 
in different spatial regions. In fact, in such cases the quantized electromagnetic field at nonzero 
temperature may not be represented in terms of (quasi)-particles [31]. As a result, the concept of a 
fluctuating field becomes not so transparent as in nonlossy dielectric media. 

In the light of this conclusion, the surface impedance approach takes on fundamental importance as 
(at present) the only self-consistent description of the Casimir effect between real metals. It does not 
need any prescription for the zero-frequency contribution to the Casimir energy and force. In all cases 
the correct expressions for the values of both reflection coefficients with two different polarizations 
at zero frequency are deduced from the general theoretical framework using the explicit form of the 
impedance function (see Sec. VI). Thus, a long discussion in the recent literature concerning the 
most adequate modification of the zero-frequency term of the Lifshitz formula [21,23-28,31-33] can 
be finalized. 

The surface impedance approach solves the puzzle with the violation of the Nernst heat theorem 
and with negative values of entropy which appears when one substitutes the Drude dielectric function 
into the usual Lifshitz formula. In the impedance approach the entropy is in all cases nonnegative 
and takes zero value at zero temperature. Thus, the general formulas given by Eqs. (41), (42), (43), 
(46), (48), and (53) lay the theoretical basis for the calculation of the thermal Casimir effect with 
respect to the needs of future precision experiments. The computations performed in Sees. V, VI are 
in good agreement with the previous results obtained by the use of the optical tabulated data and 
the plasma model. 
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The obtained results allow to remove the doubts that something is wrong with the Lifshitz formula 
[32]. In fact, the above formulas in the framework of the impedance approach coincide with the 
Lifshitz formula. The only difference is that for real metals in the frequency regions, where the 
electromagnetic oscillations initiate a real current or where the space dispersion is essential, one 
must express the reflection coefficients not in terms of the dielectric permittivity but in terms of the 
surface impedance. The usual Lifshitz formula, formulated in terms of the dielectric permittivity, 
preserves, however, major importance not only in applications to dielectrics but also in the theory 
of the non-retarded van der Waals forces between metallic surfaces. As was mentioned in Sec. II, 
the surface impedance approach is applicable with the proviso that u c < 0.1u p , i.e., the separation 
distances between the test bodies must satisfy the condition a > 5A P / (2tt) ~ \ p . In essence, the 
frequency region u > 0.1cj p is a subject of the optics of real metals near the plasma frequency [50]. At 
separations a < \ p between the test bodies the temperature effects are negligible. This is a region of 
the ultraviolet transparency where metals can be described on the same basis as dielectrics. The most 
adequate approach to the theory of the van der Waals forces at so small separations is given by the 
hydrodynamical description of an inhomogeneous electron gas [51]. This is a more general approach 
if compared with the Lifshitz theory because it does not start with a model description of a metal 
in terms of the bulk dielectric permittivity. In the local limit, however, when the spatial dispersion 
is absent, the hydrodynamical approach leads to the usual Lifshitz formula at zero temperature 
[51]. As a consequence, the usual Lifshitz formula is well adapted for the calculation of the van der 
Waals forces between real metals at separations a < \ p if e(i£) being obtained by the use of optical 
tabulated data for the complex refraction index (the extension of the avalible tabulated data into 
the region of small frequencies makes almost no effect on the computational results). 

In conclusion it may be said that the Lifshitz formula in combination with the impedance approach 
gives a solid foundation for the investigation of thermal effects onto the Casimir force. This approach 
does not lead to contradictions and can be used as the theoretical basis for the needs of future 
experiments. 
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FIG. 1. Integration paths Ci (a) and Ci (b) in the plane of complex frequency. The Matsubara frequencies are and photon 
eigenfrequencies are uj n . 
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FIG. 2. Correction factor to the Casimir energy between two Au plates at zero temperature computed by the use of the 
impedance of infrared optics (solid line) and of anomalous skin effect (dotted line) versus surface separation. 
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FIG. 3. Relative thermal correction to the Casimir free energy between two Au plates computed by the use of the impedance 
of infrared optics (solid lines) and of anomalous skin effect (dotted lines) versus surface separation at two different temperatures. 
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